Engineering non-Abelian topological memories from Abelian lattice models 
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We propose a novel fault-tolerant scheme for quantum computation based on the topological phase 
of a spin lattice model. Abelian anyons are used to simulate non-Abelian properties relevant to 
topological quantum computation. Universality is realized by using anyon braiding supplemented 
with non-topological operations, namely single spin measurements. Along with the topological 
nature of the encoding and a finite energy gap, logical states are protected from thermal excitations 
by a novel method that enforces the required degeneracies. A proof of principle scheme is proposed 
that can be realized in Josephson junctions with state of the art technology. 



Introduction:- It has been theoretically demonstrated 
that quantum computation can perform tasks that are 
virtually impossible with classical computers. Physical 
realizations are presently hindered by environmental and 
control errors that do not allow demonstration of quan- 
tum computation beyond the classical limit [l[ . Schemes 
based on topologically ordered systems promise to dras- 
tically resolve this problem by employing physical princi- 
ples as well as algorithmic constructions. Their aim is to 
encode and manipulate information in a way that is in- 
trinsically resilient to errors, thus allowing fault-tolerant 
quantum computation. 

Topological quantum computation employs anyons, 
the quasiparticles of topological models 0, 0]. Non- 
Abelian anyons have complex behavior well suited to en- 
coding and processing quantum information, but the sim- 
pler Abelian anyons are more experimentally accessible. 
Here we bridge the gap, by using well established con- 
trol techniques to engineer non-Abelian properties from 
Abelian lattice models Our method utilizes symme- 
tries of the Hamiltonian, which may be enforced by single 
spin interactions. These enable us to interpret certain 
states of Abelian anyons in terms of new quasiparticles. 
These states exhibit non-Abelian like behavior, while not 
constituting a full mapping to non-Abelian anyons. Us- 
ing these new quasiparticles, quantum information can 
be stored non- locally in the analogues of fusion spaces, 
and manipulated using single spin measurements and 
adiabatic techniques. The addition of non-topological 
operations, namely measurements of single spins in the 
underlying lattice, leads to universal quantum compu- 
tation @, la, 0, B Hi- State purification is employed to 
efficiently eliminate errors introduced by these manip- 
ulations. All employed operations, including measure- 
ments, always produce eigenstates of the Hamiltonian, 
allowing the energy gap to remain intact at all times and 
ensuring fault tolerance 10]. Though our method is gen- 
eral, we focus on the quantum double models Q, analyt- 
ically tractable topological lattice models described by 
the stabilizer formalism llj . Specifically we consider the 
D(Zp,) model, which an be realized with Josephson junc- 
tions [121 ] . An experiment is proposed to demonstrate the 



quantum memory of our scheme, as well as our method 
to enforce symmetries. 

The D(Zq) model:- The D(Zq) anyon model is defined 
on an oriented two-dimensional square lattice On 
each edge there resides a six-level spin spanned by the 
states \g) i , where g = 0, ..,5 labels an element of Zq, 
the cyclic group of six elements. The generalized Pauli 
operators for these spins are of = J2 g ez B I 9 + -0; (d I ano - 
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9~)i (9 1- We choose the orientations 
of our lattice to point upwards for each vertical link and 
right for each horizontal link. We can now define the 
following operators, 



A(v) = a*Wt of < , B(p) 



(1) 



The edges j,k 7 l,m are those sharing the plaquette p or 
vertex v. We take j to be the edge to the top of the pla- 
quette or vertex and the rest to proceed clockwise from 
this. The Hamiltonian of the model is given by, 



H = -A fe^ + EW)' 

\ V p / 



(2) 



where A > 0. A state within the ground state space, 
I gs), is then defined by the following projectors, 
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such that Pi(v) | gs) = P\{p) | gs) = | gs) for all v and p. 
All A{v) and B(p) operators mutually commute, and so 
may be described as stabilizers of a stabilizer code llj . 
The ground state space is then identified with the stabi- 
lizer space. The elementary excitations of D(Zq) are the 
anyons e g and m g for g = 1, .., 5, and correspond to viola- 
tions of the vertex and plaquette stabilizers, respectively. 
The absence of an anyon is referred to as the vacuum. 
The projectors for the anyon states are given by, 



h£Z B 

P mM = \zZ e- ivgh/3 {B(p)) h . (3) 
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When the system is in the state | ip) , a quasiparticle of 
type a at vertex v (plaquette p) is defined by P a (v) \ ip) = 
\tp) (P a (p)\ip) = \ip)). An important feature of the 
Hamiltonian @ is that it assigns the same energy, A, 
to all e g anyons, and also to all m g anyons. This sym- 
metry is not a necessary for a valid D(Zq) Hamiltonian, 
but it is necessary for the scheme we propose. Later we 
present a simple method to protect this symmetry. 

To encode quantum information in a protected way, we 
define new types of quasiparticle. We focus on the vertex 
excitations, while the plaquette excitations are similarly 
defined. Let us introduce the projectors, 

P X (v) = P e3 («), P^v) = P ei (v) + P ei (v) 
W = PeM+PeM- (4) 

In terms of the anyons of D(Zq) 1 the quasiparticle 6 is an 
e\ or e 4 and its antiparticle 6 is an e2 or eg. The quasi- 
particle A is directly identified with e^. The fact that 
the Hamiltonian assigns the same energy to each of these 
anyons allows arbitrary states of these quasiparticles to 
be eigenstates. Measuring the stabilizer code with the 
projectors in ((4]) rather than ([3]) extracts less informa- 
tion about the type of anyon present. This is equivalent 
to 'holes' in the code, where stabilizers are not enforced 
d, 0] , except that the effective holes are now carried by 
each (0, 6) pair. Our scheme uses these to store quantum 
information. 

States with (0, 6) or A pairs on vertices connected by 
a single edge, i, can be created by acting on the ground 
state with the operators: 



(5) 



respectively. The projection (1 + (of ) 3 )/2 present in Wf 
may be performed deterministically by measuring the ob- 
servable (of) 3 and applying (A(v)j to a neighboring ver- 
tex if the —1 eigenvalue is obtained. By considering the 
action of the creation operators on two edges sharing a 
vertex, one obtains the fusion rules, 



0x0=1 + A, x A = 
6x6 = 6, 6x6 = 



6 x A = 6, 
AxA=l. (6) 



In other words, if both a 6 and 6 are moved into the same 
vertex v, their state can be stabilized by either Pi(v) or 
P\(v). Note that WfWf = Wf, ensuring that the fusion 
of a A with a 6 cannot be locally distinguished from a 
single 6. This provides a non-trivial fusion space where 
information can be encoded in a non-local way. The A 
particles can be transported using chains of W x , whereas 
6 and 6 require controlled operations 0, EH such as, 



Ci= E P e g+ H(v)Pe h (v')(cr! 
g,hEZ e 
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FIG. 1: (a) Part of the lattice with edges and vertices enumer- 
ated, (b) The (0, (f>) pairs of quasiparticles in which a logical 
qubit may be stored, with the pairs highlighted. The Z basis 
may be measured by determining the fusion outcome, 1 or A 
of pairs connected by horizontal edges. The X basis may be 
accessed by single spin operations on vertical links. 



which moves the particle coherently from vertex v to the 
neighboring vertex v 1 through the edge i. Alternatively, 
our quasiparticles are well suited to be moved adiabati- 
cally using local potentials without affecting the de- 
generacy of the logical states. 

Quasiparticles x, x and /i on plaquettes can be defined 
equivalently to 6, 6 and A, respectively. The correspond- 
ing projectors P x (p), Px(p) an d Pp,{p) and the creation 
operators Wf and Wj* are obtained from ([J|, ([5]) and (JT]) 
using the substitutions A(v) — > B(p) and of — * af . The 
braiding of the quasiparticles can be determined from the 
constituent e g and m g anyons. For example, a fj, around a 
A gives the statistical phase e I7r due to their identification 
with 77i3 and e%, respectively. 

Quantum computation:- Employing the new quasipar- 
ticle states, we define a non-local encoding equivalent 
to quantum computational schemes with non-Abelian 
anyons 0, HJ ■ Consider the following operations on the 
ground state (see FigQJa) for enumeration): (i) The ap- 
plication of W^Wf; (ii) The application of Wf. Apply- 
ing operation (i) creates two (6,6) pairs at vertices V2, 
V4, and vi, V3. The fusion of either pair would result in 
the vacuum. Applying both (i) and (ii) also creates two 
(6, 6) pairs, but in this case they would both fuse to a 
A. These states can be used to encode a u-type logical 
qubit, with basis states | V ) and | l v ) associated with the 
vacuum and A fusion channels, respectively. In terms of 
spin operators acting on the ground state they can be 
written as, 



|0„) 

u 



a z 2 al[t + (alf] [t + (<7f) 3 ] | gs) , 
K)Vf oj [1 + (a*) 3 ] [1 + (o-lf] I gs) . (8) 



(7) 



No local operator (4-local such as A(v) or smaller) can 
distinguish between the two states. The fusion chan- 
nel of a (0, 0) pair is a non-local property, and can only 
be detected by observables acting on both particles of 
a pair. The logical qubit operations are given by X = 



3 



(erf) 3 or (erf) 3 , and Z = 
Consider keeping the <; 
but using ((7| to move 



A 3 ( Vl )A 3 (v 3 ) or A 3 (v 2 )A 3 (v 4 ). 
> and 4> of each pair together, 
the two pairs away from each 



other. The X operations become products of (erf) s 
along strings connecting the targeted vertices, while the 
Z operation remains the same. The two possible forms 
of both X and Z allow two possible means to measure in 
each basis. These should both give the same result, with 
any difference being a signature of errors. The fusion 
channels of x's and x's may similarly be used to encode p- 
type qubits on corresponding plaquettes pi, . . . ,pi, with 
logical states | P ) and | l p ). The X and Z rotations for 
the p-qubits are obtained using the substitutions above. 
The eigenstates of the X basis may be denoted as | ±„/ p ). 

The controlled- Z gate can be applied between a v-type 
and a p-type qubit by moving the (4>, 4>) pair on v\ , v 3 
around the (XiX) pair on p\, p 3 . This is due to the sta- 
tistical phase e l7r obtained when moving a A around a 
/i. Arbitrary single qubits rotations can be performed 
by employing suitable logical ancillary states. To do this 
we initially prepare the logical state \0 V ). The projec- 
tor H? = § [1 + cos 9 (erf) 3 +i sin 9 (erf ) 3 (crf ) 3 ] is then 
applied to the spin on edge 1. This can be performed 
probabilistically, by a measurement. In general, 11^ cre- 
ates a superposition of a pair of A's and /x's on the vertices 
and plaquettes sharing the spin on edge 1. The A's will 
fuse with the </>'s, and the [i's may be measured. If a pair 
of /i's is detected, they can be annihilated by applying 
(erf) 3 . This leaves the system in the logical ancillary state 
I a e v ) = (cos^|0„) — zsin# | l v ))/y/2. If no /x's are de- 
tected, the protocol can be repeated until successful. The 
ancillary plaquette state | ci^) can be prepared similarly. 
By utilizing these states, arbitrary single qubit rotations 
and controlled- X gates can be performed using the cir- 
cuits shown in Fig. (|2|). During the preparation stage, 
the ancilla states are stored on neighboring plaquettes or 
vertices, making them vulnerable to X errors and affect- 
ing the fault-tolerance of our scheme. However, note that 
all single qubit rotations can be constructed from e 1 ^/^ 2 
and e z (™/ 8 ^ x whose implementation requires only 



preparation of 



77/8 



tt/8 



and J +p/ v ). The former 



two may be prepared fault-tolerantly from many noisy 
copies by use of distillation 0, [TU . The latter is an eigen- 
state of X, and so intrinsically resilient to these errors. 
The controlled- AT also requires only these ancilla states, 
completing the fault-tolerant universal gate set. An al- 
ternative gate set can be performed using non-topological 
operations on the underlying spins 

Fault-tolerance:- We consider errors that do not excite 
the system, requiring a temperature low enough for topo- 
logical order to be stable IH, 17 1. The errors can then 
be considered as perturbations in the Hamiltonian, due 
to imprecise tuning of the system or coupling with the 
environment. 

The Hamiltonian ^ can be expressed as an equally 
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FIG. 2: The circuits that implement (a) e ±lBZ rotations and 



(b) 



rotations. The sign depends on the outcome of the 



measurements and can be corrected by subsequent rotations. 
In (c) the controlled- X gates are depicted on v-type qubits. 
The X measurements here may be realized by single spin 
measurements (of) 3 for «-type qubits and (af ) 3 for p-type. 



weighted sum of the stabilizers A(v) and B(p). However, 
physical systems will likely produce perturbed Hamilto- 
nians, lifting the degeneracy of the anyons and breaking 
the symmetries our scheme requires. This is a problem 
that not only affects realizations of D(Zq), but all quan- 
tum double models, Abelian or non-Abelian [3], H, Il3j |. 
Here we present a method to enforce the symmetries in 
D(Z§), but the principle applies in general. For the states 
of {4>, c6) pairs to be eigenstates of the Hamiltonian we re- 
quire both the ei and e± anyon states and the e 2 and es 
states to be degenerate. To ensure this, note that the 
operator wf involves a projection (1 + {af) 3 )/2. This 
causes the state of the pair to become an eigenstate of the 
operator (erf) 3 . Since this operation creates e% anyons, 
and so performs the mapping e\ — » e^ and e 2 — > es, it 
does not commute with any perturbation that lifts the 
required degeneracies. Application of a single spin mag- 
netic field term f?(crf) 3 then suppresses a perturbation 
of strength 5, such as SA 9 (v), by (6/B) 2 . Moving the 
quasiparticles apart means that the magnetic field is no 
longer a single (erf) 3 but a product along a chain between 
them. An equivalent method for the (x,x) pairs can be 
obtained with B(af) 3 terms. 

We now consider perturbations that do not come from 
fine tuning, such as those acting on spins forming strings 
across the lattice. If one end of a string connects with a (j> 
or <f>, it may move the quasiparticle. The simplest exam- 
ples are single spin perturbations (erf ) s , acting on a spin 
surrounding a <fi or <f>. Lowering the coupling of the ver- 
tex operators where these quasiparticles are placed from 
A — > A — a will suppress this process by an effective gap 
of size a. If string-like perturbation stretches between 
the two (</>, <j)) pairs, it can distinguish the logical states 
of the X basis and so lift the degeneracy. If the envi- 
ronment can produce string-like /c-local perturbations in 
the Hamiltonian then the ((f), </)) pairs should be moved 
k + 1 spins apart, to protect the encoded qubit. Pertur- 
bations may also act on spins that form loops. If these 
loop around a single c6 or </>, they may cause errors on 
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the results of X measurements. These will be suppressed 
by the magnetic field (erf) 3 as long as they do not loop 
around both the (j> and of a pair. Such errors become 
highly correlated, and so increasingly difficult for nature 
to produce, as these particles are moved apart. Such er- 
rors may also be dealt with by encoding each qubit on 
27V pairs rather than just 2 Q. Only perturbations act- 
ing on spins forming loops around N quasiparticles can 
then cause uncorrectable errors. 

Josephson junction realization:- We now present a 
means to experimentally demonstrate the quantum mem- 
ory of our scheme, including its resilience to perturba- 
tions breaking the required symmetries. Consider the 
Josephson junction element in Fig. Ela). A flux 2ir/6 
passes through each loop, creating six degenerate ground 
states that can be used as a six- level spin [12j . Construct- 
ing a lattice of such elements, as in Fig. EJb), imposes the 
vacuum state on all plaquettes, and gives the following 
Hamiltonian for the vertices, 



H' = -rJ2(Av)+^(v)). 



(9) 



Here 2r is the energy gap resulting from the tun- 
nelling processes within the junctions. Using a semi- 
classical approximation we find the coupling to be r rs 

^y 4 ^ /4 exp(~S*o) with So « 0.65fy^. Here Ej and 
Eq are the Josephson and charging energies, respectively. 
To realize a single qubit memory, the states (jHJ must be 
prepared. This can be done by pumping charge between 
vertices to implement the erf operations on the connect- 
ing link [18| . The degeneracy of these states must then be 
confirmed. Since the Hamiltonian is a perturbed version 
of (J2J), it does not assign equal energies to the e g anyons. 
Hence, the terms i?(erf ) 3 and -B(erf ) 3 must be applied to 
each pair to ensure the degeneracies. These can be sim- 
ply implemented by passing a flux of 27r/3 through the 
elements on links 2 and 3, rather than 27r/6. The ex- 
perimental verification of such memories, even without 
the application of any quantum gates, would form a ma- 
jor breakthrough in the realization of anyonic quantum 
memories. 

Conclusions:- We have demonstrated non-Abclian like 
encoding of quantum information using the more exper- 
imentally accessible Abelian anyon models. Since the 
scheme requires a fine tuned Hamiltonian, we have also 
introduced a method to enforce symmetries that can be 
applied to other models, both Abelian and non- Abelian. 
Further, we have proposed an experiment to demonstrate 
the quantum memory and enforcing of symmetries with 
cutting edge technology, ft would be interesting to study 
whether the enhanced fault-tolerance provided matches 
that of non-Abclian schemes requiring non-topological 
operations for universality la ]. 
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FIG. 3: (a) Twelve Josephson junctions arranged in five loops. 
This element has six degenerate ground states when a flux 
27r/6 is passed through each loop, and thus provides the re- 
quired six-level spin, (b) A plaquette of the D(Z&) model, 
equivalent to that shown in Fig[]Jb) , realized with Josephson 
junction elements. 
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